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Abstract 

We obtain exact expressions for the two-time autocorrelation and 
response functions of the d-dimensional linear Glauber model. Al- 
though this linear model does not obey detailed balance in dimensions 
d > 2, we show that the usual form of the fluctuation-dissipation ra- 
tio still holds in the stationary regime. In the transient regime, we 
show the occurence of aging, with a special limit of the fluctuation- 
dissipation ratio, = 1/2, for a quench at the critical point. 
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1 Introduction 

It is widely recognized that the ferromagnetic Ising chain with first-neighbor 
interactions and Glauber dynamics is one of the simplest, exactly soluble, 
stochastic dynamical systems jl]. At finite temperatures, in the stationary 
regime, the two-time spin autocorrelation, C(t,t r ), and the associated re- 
sponse function, R (t, t'), of this Glauber chain are time-translationally invari- 
ant, and duly related by the usual expression of the fluctuation-dissipation 
theorem [2]. There is also an aging regime, with violation of the usual form 
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of the fluctuation-dissipation theorem. At the critical point, at zero temper- 
ature, for large values of the observation time t, the fluctuation-dissipation 
ratio of the Glauber chain assumes the non trivial limoting value = 1/2. 

We now revisit a linearized version of the Glauber model, proposed by 
one of us a few years ago 3J. As in the original Glauber model on a d- 
dimensional hypercubic lattice, we still consider a one-flip stochastic process. 
Each site r = 1, N of the lattice is associated with a spin variable o r = ±1. 
However, the time evolution is now governed by a linear spin-flip ratio, 



where AG [0, 1] is a parameter, the sum is over the 2d nearest neighbors of site 
r, the time scale is set by the parameter a, and a stands for a configuration 
of spin variables, a = {<J r }. The evolution of the probability P(a,t) of the 
spin configuration o at time t is given by the master equation 



where a r is defined as the configuration a with a r replaced by — a r . From 
these equations, it is straightforward to calculate exact analytical expression 
for the site magnetization, m r (t) = (a r (t)), and the pair correlation function, 
Qr,r'{t) = (&r{t)o'r'{t))- In contrast to the original Glauber model, defined by 
a non linear transition rate, with exact solutions restricted to d — 1, we 
can now write expressions for m r (t) and q r ,r'{t) in all dimensions d. For 
< A < 1, the linearized model displays a disordered (paramagnetic) phase, 
with exponentially decaying pair correlations; at the critical point A = 1, 
correlations decay algebraically [3]. 

It is important to distinguish between this linear Glauber model and the 
so-called Glauber dynamics, associated with a nonlinear transition rate, and 
which is used to simulate the Ising model. The linear Glauber model may 
be regarded as a voter model with noise |3]. It displays just one phase, for 
all dimensions, as long as noise is finite. In the absence of noise (A = 1) 
it becomes critical. In one dimension there is a close identification between 
these two versions of the Glauber model. For d > 2, however, in contrast to 
the original model, the analytically solvable linear Glauber model is micro- 
scopically irreversible (in other words, although having a stationary state, it 
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dt 



P(a, t) = J2 KK)PK, t) - w r (a)P(a, t)} , 



(2) 
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does not obey detailed balance, and cannot be associated with a Hamilto- 
nian). The fluctuation-dissipation theorem is usually conceived for systems 
that do obey detailed balance E] ■ It is then reasonable to ask some ques- 
tions, including the validity of the fluctuation-dissipation theorem and the 
presence of an aging regime, about the dynamical behavior of systems that 
do not obey detailed balance. One of the purposes of this article is to carry 
out a thorough analytical investigation of a particular system, as the non 
linear Glauber model in d > 2 dimensions, which belongs to the large class 
of microscopically irreversible models [3J EJ IHl E] • 

The dynamical calculations of interest in this investigation are performed 
in the presence of a (small) perturbation. In the treatment of stochastic 
models it is natural to introduce the modified one-spin-flip rate, 

Wr (a) = w° r (a)e- hra % (3) 

where w° (a) is the unperturbed flipping rate associated with the rth spin, 
and h r is the time-dependent disturbance coupled to the dynamic variable 
<j r . Alternatively, as the calculations are restricted to small perturbations, 
we can write 

w r (a) = w°{a){l - h r a r ). (4) 

If the model is microscopically reversible, that is, if the unperturbed transi- 
tion rate w° obeys detailed balance, there is a model Hamiltonian 7i , and 
it is straightforward to show that the perturbed transition rate w r , given by 
Eq. (J3J), also obeys detailed balance. In this reversible case, the model is 
described by the Hamiltonian H = H a — J2 r H r a r , where H r = h r / j3, and is 
/3 proportional to the inverse of the temperature. This form of perturbation 
is then suitable for reversible models, with a disturbance h r proportional 
to the external field. Assuming this expression for the flipping rate, the 
fluctuation-dissipation theorem is given by 

R(t,t') = -^ f C(t,t% (5) 



where 

5m r (t) 



N ^ Sh r (t') 



(6) 
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is the response function, m r (t) is the average of a r at the observation time t, 
and 
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is the autocorrelation function of a r between the observation time t and the 
waiting time t' (with t > t'). For systems obeying detailed balance, we have 
h r = /3H r , and relation (JBJ) reduces to the usual form of the fluctuation- 
dissipation theorem. Another version of the fluctuation-dissipation theorem 
relates the susceptibility, associated with a given dynamical variable M, such 
as the total magnetization, and its variance, 

J-(M) = (M 2 )-(M) 2 , (8) 

where h is a static, time-independent disturbance, introduced by the pre- 
scription of Eq. (JHJ). 

In this work, we show that both forms of the fluctuation-dissipation re- 
lation, (jSJ and (|8|). are valid for the linear Glauber model in the stationary 
regime. In the transient regime, where aging behavior takes place, these 
relations are no longer obeyed. It is then appropriate to define |1U1 ITT] a 
fluctuation-dissipation ratio, 

In the linear Glauber model, for all values of the dimension d, we show that 
X(t, t') — ► 1 in the limit if — > oo, except at the critical point, A = 1, in which 
case X(oo,f) -> 1/2. 

The layout of this paper is as follows. Some results for the linear Glauber 
model, including a discussion of the lack of detailed balance for d > 2, and 
calculations of the site magnetization and the spatial two-body correlations, 
are reviewed is Section 2. These one-time functions play a major role in 
the calculations of the two-time functions, as the autocorrelation and the 
response functions, which are obtained in Section 3. In this Section, we 
also consider the stationary limit and make a number of comments on the 
non- stationary regime. Section 4 contains some conclusions. 



2 The linear Glauber model 

We have already mentioned that the linear Glauber model has been intro- 
duced by one of us [3| as an extension of the original Glauber model [I], and 
may be regarded as a voter model with noise. This linear Glauber model is 
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defined by the (linear) one-spin- flip rate given by Eq. (JTJ, which should be 
inserted into the master equation (J5J). 

^From equations (JTJ) and (J2J), it is not difficult to obtain the evolution 
equations for the site magnetization and the spatial pair correlation, 

~dt mr ^ = ~ mr ^ + 2d 5^ mr + 5 ^ ( 10 ) 

8 

and ^ 

--rAr,r'{t) = -2q rr ,(t) + — V" [g r ,r'+<$(*) + g r ',r+a(*)] , (11) 
a at la *—r 



for r / r'. If r and r' are nearest-neighbor sites, the right-hand side of Eq. 

contains terms like q r , r (t) and q r 'y(t) which should be set equal to 1. 
The possibility of obtaining these exact expressions, and of performing the 
exact calculations that are going to be reported in this article, is one of the 
most relevant features of the linear Glauber model. As shown by Oliveira [3], 
for < A < 1 this model displays a disordered (paramagnetic) phase with 
exponentially decaying correlations. For A — > 1, it becomes critical, with 
algebraically decaying correlations at A = 1. 

The linear Glauber model in one dimension has a reversible dynamics. In 
one dimension, the probability of occurrence of any sequence of states and 
the probability of the associated reverse sequence of states are the same. In 
dimensions larger than one, this is no longer valid. Consider, for instance, 
the four states shown in Fig. [TJ on a square lattice. Suppose that the system 
follows the sequence of states A, B, C, D, and returns to the initial state 
A. If the interval At between two successive states is small, then, according 
to the spin-flip rate given by Eq. Q, the probability of occurrence of the 
sequence A— >B— >C— >D— >A is given by 

P(A^B^C^D^A) = P{A\D)P(D\C)P{C\B)P{B\A)P{A) 

= ^(l-^ 2 (l + A)(aAt) 4 P(A). 

On the other hand, the reversed sequence, A^D^C-^B^A, has the prob- 
ability 

P(A —>D—>C—>B^A) = P{A\B)P{B\C)P{C\D)P(D\A)P{A) 

= ^(l + ^ 2 (l-A)(«At) 4 P(A). 
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Figure 1: A possible irreversible sequence. 



These two probabilities are distinct, so that the linear Glauber model on a 
square lattice is indeed irreversible, the only exception being the trivial case 
A = 0. A generalization of this result to larger dimensions is easily found 
by filling the sites created by the introduction of more dimensions with "+ 
spins" . Hence, the detailed balance cannot be valid, and the stationary state 
is a priori not known. The connection of the transition rates with a Gibbs 
measure, as it has been possible in the one-dimensional case, is now forbidden. 

2.1 Site magnetization 

We now introduce the Fourier transform of m r (t), 



Jo 

Using these transforms, the differential equation (jTUjl is reduced to the alge- 
braic form 





r 



and the Laplace transform of rrik(t), 




(13) 




(14) 



where 




(15) 
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and rh° k is the Fourier-Laplace transform of m r (0). 
The inverse Laplace transformation leads to 

m k (t) = e- a ^m fc (0). (16) 

Inverse Fourier transforming, we then have 

m r(t) = / ] r r - r /(t)m T ./(0), (17) 



where 

rr(t) = / e'-'-^'gj. (18) 

We can also obtain the site magnetization with reference to an initial time 
t = t' instead of t — 0. In this case, we just write 

m r (t) = ^T r _ r ,(t-t')m rl (t'), (19) 

r' 

where m r (t') is the site magnetization at time t' . 
2.2 Pair correlation 

If we look for translationally invariant solutions of Eq. (fTTJ) . the spatial 
correlation between sites r and r' will be a function of distance, (cr r (t)cr r i (t)) = 
q r - r '{t) = q r i_ r (t). We then write Eq. (fTTj) as 

8 

for r ^ 0, with qo(t) = 1 whenever it appears on the right-hand side. 

Using a method introduced by Oliveira [3], let us write an equation for 
r = 0, 

8 

where b(t) is chosen to ensure that qo(t) = 1. Actually, b(t) is defined by 

& (*) = 2- -£>(*). ( 22 ) 

5 
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Consequently, Eqs. (J2*U|) and (|2*T]) can be written as 



for all values of r. 

Eq. (J23j) can be examined from two points of view. On the one hand, 
it is a first-order differential equations in time. On the other hand, the 
summation over nearest neighbors resembles a discrete lattice Laplacian, and 
from this point of view it is a discrete second-order difference equation that 
can be solved by the use of a Green function. We then introduce the Laplace 
transform of q r (t), 

POO 

q r (s) = / dte- st q r (t), (24) 
Jo 

so that Eq. (j2*B"j) may be written as 

- [-ml + sq r (s)] = -2q r (s) + 3 q r +s(s) + b(s)S r>0 , (25) 

o 

where b(s) is chosen such that q (t) = 1, which is equivalent to take g (s) = 
1/s. For a random initial condition, which corresponds to a quench from 
infinite temperature, it is appropriate to take q r (0) = ml (1 — <5 ri0 ) + ^r,o- 

To solve equation (|23|) . we introduce the lattice Green function associated 
with a d- dimensional hypercubic lattice, 

/ ikr Jdr, 
J7TT1 w . ( 26 ) 

where the integration is over the first Brillouin zone, and 

A d 

/(A0 = l--5>os%. (27) 
3=1 

This Green function satisfies the equation 

-sG r (s) = -2G r (s) + ^ d r+s(s) + S rfi . (28) 

o 
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In terms of this lattice Green function, the solution of Eq. (|15|) can be 
written as 



rrir 



q r {s) 



+ b{s)G r {s) 



s + 2a(l - A) 

where b(s) should be chosen so that qo(s) = 1/s. This leads to 



b(s) 



G (s) 



_s s + 2a(l-A) 



(29) 



(30) 



from which it follows [3] the solution 



Qr(s) 



TTln 



s + 2a{\-\) 



1 



G r ( S ) 



+ 



G r (s) 

sG (s) 



(31) 



For a completely random initial condition, m = 0, we have m{t) = 0, 

1 



and 



b(s) 



Qr{S) 



sG {s) 

G r (s) 
sGn(s) 



(32) 



(33) 



3 Two - time response and autocorrelation 
functions 



The calculation of the two-time response function 



N ^ 5hJt>) 



(34) 



requires the application of a small perturbation, which is introduced accord- 
ing to the prescription P|). and from which one measures the response of the 
system. As pointed out in the Introduction, we assume a perturbed spin-flip 
rate, given by 



t \ a 
w r (a) = - 



A 



[1 — h r {t)a r 



(35) 
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where h r (t) is a time-dependent disturbance coupled to the dynamic variable 
a r (t). We then write the equation of motion for the site magnetization, 



-^-m r (t) = -m r (t) + ^^m r+& (t) + h r (t) 



A 



(36) 



which can also be written as 

-m r (t) = -m r {t) + ^- y^m r+& {t) + ^h r (t)b(t), 



1 d 

a dt 



(37) 



where is given by Eq. 

Using the same procedures adopted to solve Eq. (JTUJ), and taking into 
account that the last term in Eq. (|37|) is a known function of time, we have 



m 



.(*) = Vr r _ r /(£)m r /(0) + / ^r-r'(t-t')h rl (t')b(t')dt'. (38) 

2 „, ^0 



^From this expression, we calculate 

Sm r (t) a 



5h r ,(t') 2 
which leads to the response function 



r r _ r ,(t-t>(t'), 



R(t,t') = -r (t-t')b(t'). 



(39) 



(40) 



The correlation (a r (t)a r (t')) of a spin at a given site r, at time t', with 
the same spin at a later time t (t > t') is formally written as 

(a r (t)a r (t')) =Y,Y, a r(t)PM<7 , ,t'K(t')P(<r , ,t'), (41) 



where P(a,t\a',t') is the conditional probability of finding the configuration 
a at time t given the configuration a' at an earlier time t' . Noting that the 
site magnetization, m r (t) = (a r (t)), with the initial condition m r (t') = cr' r (t'), 
may be written as 

<7 r P(v, t\cr', t') = m r (t), (42) 
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and using Eq. ()19j). we have 

r (t) = ^ y r r _ r /(t — t')a r '(t'), 



m r 



which can be inserted into Eq. (J41)) to give 

(a r (t)a r (t')) = J2J2 T r-r>( t - t >'At'KP(v , ,t'), 

a' r' 

which finally leads to 

C(t ) t') = Y,?r(t-t')q r (t'). 

r 

3.1 Stationary regime 

In the stationary regime the waiting time t 1 and the observation time 
without limits, but the difference £ — t' is fixed. To be more precise, t' 
with t > t' and r = t — t' fixed. 

^From the Laplace transform final value theorem, we have 

q r (oo) = lim q r (t) = limsg r (s) = 

s ~*° G (0) 

so that 

C (t,t') = C(r) = Vr r (r)^, 
V Go(0) 

which can be written as 

C(r) = -J— / -L^- ^ — 



G (0)7 V\k) (27T)-" 
On the other hand, taking into account that 

6(00) = lim = limsS(s) = — , 

the response function ()4()j) can be written as 

a 1 
i*(t,0 = #(T) = irro(r 



2 'G (0) 
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Using the definition of T r (t) , given by (|TH|) , we have 



R ^^l'-" mT Wf- (51) 

Both quantities, C(r) and R(t), are time-translationally invariant (func- 
tions of t only), as it should be anticipated in a stationary regime. More- 
over, the fluctuation-dissipation theorem is trivially satisfied, with R(t) = 
—dC(r)/dr. For large time differences, and A 7^ 1, it is easy to see that both 
the autocorrelation and the response functions decay exponentially, accord- 
ing to exp [—a (1 — A) r], with the equilibration time r eq — l/[ot (1 — A)]. 

3.2 Global fluctuation-dissipation relation 

The fluctuation-dissipation theorem can also be written in terms of global 
variables, as the magnetization and the corresponding answer to a static 
perturbation. Let M denote the magnetization, 

r 

and let us consider a static homogeneous disturbance h, defined by Eq. (jSJ). 
Then 

dM(t) 

-or = (53) 

where x{t) is a variance, 

*(*) = ^EE^W) " M*)X<v(f)>l- (54) 

r r' 

Due to the translational invariance of the lattice, M(t) = (<7o(i)}. This 
quantity is the solution of Eq. (pT7|) . which is now given by 

i^M(t) = -(l-A)M(t) + i^(t), (55) 
a at 2 

with the stationary solution 

M = h-^- = h . (56) 

2(1 -A) 2(1-A)G (0) 
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For A < 1, the magnetization vanishes as h — > 0, so that the variance is 
given by 

*(*) = $>(*). (57) 

r 

^From Eq. (J23j) it follows that 

= "2(1 -A) X (t) + &(«), (58) 

with the stationary solution 

= Hop) = 1 

2(1 -A) 2(1-A)G Q (0)' 1 ' 

^From equations (JoT)j) and ()59p. it is seen that relation (|53"j) is clearly satisfied 
in the stationary regime. 

3.3 Aging regime 

The temporal behavior of the autocorrelation and the response functions, as 
calculated in this section, already suggests the existence of an aging regime. 
Equations and (|45p . for R(t,t') and C(t,t'), are valid for all values of 
t' and t, with t — t' = r > 0. The dependence of these functions on both t 
and and not on r only, leads to the existence of aging. In this regime, the 
role of the spatial correlations is still crucial, since they are responsible for 
the realization of the aging scenario. In the stationary regime, it should be 
noted that q r becomes a time-independent quantity in the t' — > oo limit only, 
and this is the reason of the dependence of the two-time functions on r only. 

In the transient regime, we do not expect the validity of the usual form of 
the fluctuation-dissipation relation given by Eq. (jSJ). It has been convenient 
[TU1 ITT| El to characterize the distance to the stationary regime by the 
fluctuation-dissipation ratio, 

A particular interesting quantity is the limit 

= lim X(oo,t'), (61) 

f— +oo 
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where 



X(oo,t') = lim X(t,t'). 



(62) 



^From the results of this Section, it is easy to write 




" 2b{t')- 2(1 -A)x(t') 



(63) 



where x(0 i s the sum given by (|H7j) . At the critical point, A = 1, it follows 
that X(oo,tf) = 1/2 for any time t', so that Xoo = 1/2. In the disordered 
phase, A 1, Eqs. ((SHI) and (|55Jl may be used to conclude that = 1. 

4 Conclusions 

We have reported a number of exact calculations for the dynamical behavior 
of a d- dimensional linearized version of the stochastic Glauber model. In one 
dimension, both the linear and the original model are essentially equivalent. 
For d > 2, however, the rates of transition of the linear Glauber model 
do not obey the conditions of detailed balance. This linear model can be 
regarded as a voter model with noise; it displays just one stable phase as 
long as noise is finite, and becomes critical in the absence of noise (A = 1). 
Since the dynamical properties are usually conceived for systems that do 
obey detailed balance, we decided that it was appropriate to carry out some 
explicit calculations for a microscopically irreversible model. 

We have obtained expressions for the two-time autocorrelation, C (t,t'), 
and response functions, R(t,t'), which depend on both observation time t 
and waiting time t' < t. In the stationary, infinite time limit, t' — > oo, in 
the presence of noise (0 < A < 1), the spatial correlations are independent of 
time, and the two-time functions become translationally invariant (depending 
on the difference r = t—t' only). The usual form of the fluctuation-dissipation 
theorem, X(t,t') = R(t,t')/dC(t,t')/dt' = 1, is trivially observed in this 
regime. 

In the scaling regime (t — > oo), for < A < 1, we obtain a non trivial 
fluctuation-dissipation ratio, X(oo,t') ^ 1. At the critical point, A = 1, we 
have X(oo,t') = 1/2, which further indicates that the dynamical behavior of 
the c?-dimensional microscopically irreversible linear Glauber model is very 
similar to the behavior of its one- dimensional reversible counterpart. 
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